Quantum computation via measurements on the low temperature state of a 

many-body system 
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We consider measurement-based quantum computation using the state of a spin-lattice system in 
equilibrium with a thermal bath and free to evolve under its own Hamiltonian. Any single qubit 
measurements disturb the system from equilibrium and, with adaptive measurements performed at 
a finite rate, the resulting dynamics reduces the fidelity of the computation. We show that it is 
possible to describe the loss in fidelity by a single quantum operation on the encoded quantum state 
that is independent of the measurement history. To achieve this simple description, we choose a 
particular form of spin-boson coupling to describe the interaction with the environment, and perform 
measurements periodically at a natural rate determined by the energy gap of the system. We found 
that an optimal cooling exists, which is a trade-off between keeping the system cool enough that 
the resource state remains close to the ground state, but also isolated enough that the cooling does 
not strongly interfere with the dynamics of the computation. For a sufficiently low temperature we 
obtain a fault-tolerant threshold for the couplings to the environment. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

The "one-way" model for quantum computation, 
which requires only local adaptive measurements of in- 
dividual qubits prepared in a fixed multi-qubit resource 
state, provides a new approach for assessing the physi- 
cal requirements for universal quantum computing. The 
cluster state on a two-dimensional square lattice is 
the canonical example of a resource state that allows 
for universal measurement-based quantum computation 
(MBQC) P, Much research on MBQC focusses on 
properties of the resource state itself, and in particu- 
lar how such a state could be prepared dynamically via, 
say, local controUcd-Z operations in a variety of sys- 
tems for which the dynamics of the individual qubits 
can be uncoupled, such as an optical lattice [1] or sin- 
gle photons [1, [Hi- In contrast, recent new theoretical 
results in MBQC have been obtained by viewing the 
resource state as the ground or low-temperature ther- 
mal state of a strongly-coupled quantum many-body sys- 
tem H 0, i, i, lia, nil. This perspective allows us 
to use some powerful tools and techniques of quantum 
many-body theory, for example, to determine what type 
of systems permit universal MBQC [13, 0, ES] and for 
those that do, how robust the system is in its universal- 
ity iiiP- 

One could take this perspective of ground states serv- 
ing as computational resources as a physical realisation, 
and thus obtain a new mechanism for creating cluster 
states or other such resource states. That is, if a quan- 
tum many-body system could be engineered such that it 
possesses the cluster state as its unique ground state , 
and if the system is sufficiently gapped, then MBQC can 



be performed by cooling the system down to a sufficiently 
low temperature and then performing a sequence of adap- 
tive measurements on this thermal state. However, by 
treating the resource state for MBQC as the equilibrium 
state of a dynamical system, any measurements that we 
perform will necessarily disturb it from its thermal state. 
As measurements must be adaptive, if they are separated 
by finite time intervals we are faced with both errors pro- 
duced by the evolution under the system's Hamiltonian 
and also the cooling interaction with the environment. 
These two sources of dynamical error, together with ther- 
mal errors, act to reduce the output fidelity of any MBQC 
scheme that we might wish to perform. 

In this paper, we investigate a regular lattice of qubits 
for which the free Hamiltonian has the cluster state as its 
ground state. The system is first prepared by cooling via 
a simple and convenient choice of coupling to a bosonic 
bath in a thermal state, and we assume that the coupling 
to the bath is present throughout the computation. This 
situation is relevant to an experiment in which a strongly 
coupled system is first prepared in a useful initial resource 
state using a refrigerator, and which cannot easily be 
subsequently decoupled from the refrigerator before the 
MBQC commences. Alternatively, in the context of a 
laser-cooled atomic system, it may be inconvenient or 
undesirable to turn off the cooling interaction before the 
MBQC commences. 

We explicitly determine how MBQC proceeds on this 
system's thermal state as it is perturbed by measure- 
ments, with free evolution and cooling interaction ongo- 
ing between measurements. In particular, for the lattice 
of spins in the presence of a spin-boson coupling to a 
thermal bath that acts to restore the system to the pure 
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cluster state, we show that the free Hamihonian for the 
spin lattice state induces coherent oscillations that de- 
termine a natural measurement rate. Importantly, the 
effect of the bath may be conveniently described by a sin- 
gle quantum operation that acts on the encoded quantum 
information within the MBQC computation and which is 
independent of the particular measurement history. With 
this result we demonstrate that MBQC on such a dynam- 
ical thermal state is fault-tolerant for a sufficiently low 
temperatures and for couplings to the bath below a given 
threshold. 



II. MBQC WITH DYNAMICS 

The cluster state on a lattice C oi N qubits may be 
defined in the stabilizer formalism [isj as the unique 
eigenstate of each of the mutually commuting stabilizers 
Ki = Xi (Sij^i Zj with eigenvalue one 0] ■ Here i labels a 
particular site in the lattice and j ~ i signifies that j is a 
neighbouring site of i. The stabilizer description allows 
us to define the cluster state as the ground-state of the 
Hamiltonian 

iec 

with an energy gap A. A state in the fc"" excited level 
of this Hamiltonian is obtained by performing Z errors 
at k distinct sites, and implies that the energy level is 
(^)-fold degenerate. 

A useful alternative description of the cluster state, 
which we shall make use of shortly, is in terms of the ac- 
tion of entangling unitaries between neighbouring qubits 
on the lattice. The qubits at all sites are first initialized 
in the state |-|-), i.e., stabilized by the set of operators Xi, 
and then for every bond in the lattice the controlled-^ 
unitary exp(z7r|l)(l| (g) |1)(1|) is performed between the 
two end qubits. We denote the product of all controUcd- 
Z operations on each bond simply as U, and the link 
between the two descriptions is provided by the relation 

UK,W = X^. 

MBQC on the ground state of a spin lattice model gov- 
erned by this Hamiltonian involves an adaptive measure- 
ment procedure, in which qubits are measured sequen- 
tially in different bases until the desired output state is 
produced, up to Pauli operator corrections, on the re- 
maining unmeasured qubits. The Z measurements play 
a special role and are used to remove individual qubits 
from the cluster state, while sequences of single qubit 
measurements in the X-Y plane are used to implement 
unitary gates on the encoded quantum information. For 
the perfect cluster state, the inputs can be taken, with- 
out loss of generality, to be |+) on each of the qubits to 
be measured first. 



A. Measurements and free evolution of the lattice 

As some of these measurements arc adaptive (i.e., the 
choice of measurement bases are conditional on prior 
measurement outcomes) they must necessarily be per- 
formed at different times. Measurements disturb the sys- 
tem out of its ground state, and between measurements 
this disturbed state will evolve under the Hamiltonian 
P]). An important property of this Hamiltonian is that 
it is dispersionlcss, and so any localized excitations will 
remain local and will not propagate across the lattice. 
For example, if a Z measurement is performed at site 
s the system is projected into an equal superposition of 
the ground-state U and the state with a single Z 
error on site s, U{(^i^s\+)i ® K)s)- In the case of a sin- 
gle X measurement at the site s, the system is projected 
into a superposition of the ground-state U{^i\+)i) and 
the state with Z errors on all of the neighbouring sites 
of s, U{iE)i^s\+)i ^ir^s This local disturbance re- 

mains local under evolution; however, because the post- 
measurement state is no longer an energy eigenstate, this 
evolution must be accounted for when we perform sub- 
sequent measurements on neighbouring qubits. 

If all of the measurements on the system can be per- 
formed on a timescale much less than that of the sys- 
tem's evolution, one may be able to treat the effect of 
short-time evolution as a small perturbation of the clus- 
ter state. However, one could alternatively make use of a 
natural timescale of this system. The equal-spaced spec- 
trum of the Hamiltonian ([T]), with spacing A, ensures 
that evolution is periodic with period t = 2tt/A. If 
(instantaneous) measurements are made at time inter- 
vals which are multiples of this period, the evolution of 
the system under the Hamiltonian can be ignored. In 
essence, the energy gap A of the Hamiltonian provides 
a natural "clock speed" for the quantum computation. 
Given that the gap in the system will determine the tem- 
perature to which the lattice must be cooled in order to 
approach the ground state, it will be desirable to engi- 
neer systems in which this gap is as large as possible; with 
this in mind, performing very fast measurements (i.e., at 
a frequency S> A) may not be possible, and performing 
measurements at this clock speed (or integer fractions 
thereof) is a much less stringent requirement. 

B. Spin-boson model 

The quantum many-body system with Hamiltonian ([T]) 
is gapped, and so we can prepare a cluster state (or a 
close approximation to it) by cooling the system down 
to near its ground state though coupling to a thermal 
bath. (This cooling can be done efficiently because of 
the simplicity of the Hamiltonian (l6j.) Performing mea- 
surements on the ground state yields excited states that 
are no longer in equilibrium with the bath and so, if the 
cooling interaction is present, any measurement scheme 
that we may perform on the cluster state must proceed 
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sufficiently quickly to avoid a return to equilibrium. How- 
ever, we have already argued that the free Hamiltonian 
He will require measurements to be close to the intervals 
27r/A, and this clock speed provides a lower bound on the 
overall duration of the computation. We now consider 
the effect of a finite measurement rate in the presence of 
such cooling. 

To model the effects of cooling, we consider a system 
consisting of a bath of bosons held at a low temperature 
and coupled to each site qubit via a spin-boson interac- 
tion, which takes the form Hj = J2i j ■^ij ^ilij ■ Here, 

Qij = o.ij + alj is the displacement operator for the j"* 
mode at site i and Ay are coupling constants. The full 
Hamiltonian for the system of qubits and bosons is then 



Hi + Hb , 



(2) 



where Hi, = '^^j uJijaljaij is the free Hamiltonian for the 
bath. 

We note that our results depend on this choice of Z axis 
to describe the coupling to the bath. In practice we may 
not have full control over this coupling, but in many sys- 
tems (e.g., in trapped atoms), to a good approximation 
the environment couples only to a single spin component 
of the qubit degree of freedom. In such a situation, one 
may take the coupling to the bath as defining the Z axis. 
This assumes we have full control over the cluster Hamil- 
tonian, and so may adjust it so as to coincide the Z axis 
for the cluster state with the axis defined by the cooling 
interaction. We leave as open the question of how MBQC 
can proceed with a more general coupling. 

Because the interaction Hamiltonian commutes with 
the set of controUed-Z unitaries applied to every neigh- 
bouring pair of qubits. we can map this system using the 
unitary [/ to a dual system of uncoupled qubits, with the 
same interaction Hj to the thermal bath. We shall con- 
sider the master equation for this dual system with total 
Hamiltonian 



if. 



Hi + Hb 



(3) 



In general, we denote operators in the dual model with 
a overline bar. For example, single-qubit measurements 
on the original system, given by projectors Pk, are now 
described in this dual model by multi-qubit projectors 

'K=uPkW. 

A standard derivation results in a master equation [l7j 
for the lattice subsystem given by 

^(*) = E (4 ^(*)] + [I +>. (- 1] m 



+ p,v[\~y,{+\]pit)), (4) 

where the action of the superoperator 2? [A] on the state 
p IS given by V[A]p = ApA'< ~ ^{A'<A,p}. The constants 
ai and Pi are parameters that depend on the couplings to 



the bath. Ay, and the temperature of the bath. T. They 
are given explicitly as 

a, = 27r^A?/K-A)(l-^n(A)) 

j 

A = 27r^A|.(5K - A)7i(A) 



i{E) = (e 



J 



(5) 



We make the simplifying assumption that the couplings a 
and P do not vary from site to site and for later reference 
we may relate the temperature of the bath to the coupling 
parameters through the equation 



kT = A/(log(a//3)). 



(6) 



Within the dual picture, a qubit initially in the state 
Po will evolve in time t under a completely-positive (CP) 
map St to the state = £t {7>o)- ^ Kraus decomposition 
for this CP map £{p) — J2i ^i.tpM-^, is given by 



Ml. = 



a 



-iAi 



a + P 



-){- 



-(a+/3)t/2| 



-)(-|) 



M- 



2,t 



P 



a + P 



(e+^'^V)(-| + e-("+^^*/'|+)(+|) 



Ms,* = y^a(l-e-("+«*)/(a + /3)|+)(-| 

M4,t = y'/3(l-e-("+«*)/(a + /3)|-)(+| . (7) 

This evolution takes any single qubit state p asymptoti- 
cally in time towards an equilibrium state 



a + P' 



-){- 



P 



a + P 



1 



1 + g-A/fcTl 



-)(- 



l-X-l 

1 + e-^/^T I 



-X-l- (8) 



Thermal equilibrium for the full lattice is reached with a 
rate governed by the couplings a and p. 



C. Example: Arbitrary X-rotation 

To illustrate the effect of dynamics on MBQC we con- 
sider performing a simple single-qubit gate using MBQC 
on a one-dimensional lattice. More general gates will be- 
have similarly, as we shall show in Sec. IIIII 

Consider performing an arbitrary X-rotation gate, 
i.e., a rotation X{6) = exp[— of a single qubit 
about the X-axis by angle 9. The smallest clus- 
ter state that can realise such a gate is the thrce- 
qubit cluster state on a line. The ideal gate pro- 
ceeds as follows for a non-dynamical cluster state. 
The qubits are initially prepared in the state \ipin)i ® 
+}2 ® l+)3- The state is then entangled with the uni- 
tary U. Qubit 1 is measured in the basis {|+),|— )}, 



4 



with measurement result si € {0,1}- Based on this 
measurement result, qubit 2 is measured in the ba- 
sis {exp(-ifZ)|+),cxp(-i§Z)|-)}, where = {-lY'O, 
with measurement result S2- For the static case it is 
simple to show that, subsequent to these measurements, 
qubit 3 is left in the state exp[—i^Xs]Z^^X^'^\^)s = 
^■3^X3^ exp[-i^X3]\ip)3. That is, the initial state lip) has 
been rotated by the gate exp[— up to Pauli operator 
corrections Z^^X^^. 

For a dynamical thrce-qubit cluster state that evolves 
according to the Hamiltonian ([2]), the timing of the two 
projective measurements becomes important for the gate 
to succeed with high fidelity. First, if the initial state is 
left to interact with the bath, the system would eventu- 
ally evolve to the equilibrium state and the input state 
An = |'/'in)(V'in| wouH bc crascd. For temperature T > 
we assume that the initial state is U{p.^®'f^'^)U\ where 
Pg is given by Eq. ([5]), and that the system evolves for 
a time to until the projective measurement Pi on qubit 
1; the measured state then evolves until time to + t a.t 
which point the second measurement P2 is performed. 
The output state is thus given by 

p = UP2£t(Pl£to (Pin ® pf')P\)Plu^ , (9) 

where we described the evolution in our dual model, re- 
lated to our system by the unitary operation U, with evo- 
lution at time U obtained from (0]), and Pk = UPkU^ . 
The evolution up to time to is given by 

fto :ft„®pf -fto(pJ®pf . (10) 

For convenience, we define p^^^ — ^to(Pi„)i which can be 
expressed in Bloch vector form as = i(^ + '^to ''^) with 
rto = {xto,yto,Zto)- 

The second stage of the evolution is different due to 
the projective measurement on the first qubit. A direct 
calculation of ([9]) followed by tracing out qubits 1 and 2 
yields the final state of qubit 3 

Mto,t) = Zl^X^^e~^-^^^p,,,e^^^^-Xl^Z-;^ , (11) 

where ^ ^[I + ^out(io, t) ■ ^] and 

n>ut(to,i) = (a;out, yout,-Zout) 

= a;toe"^("+'^)*cosAt 
y^^^ = (yj^ COS At - zta siu At)e~("+'^^* cos At 
Zout = (zto cos At + yt„ sin At)e"5("+/5)t , (12) 

The decohcrencc due to the evolution under the coupling 
to the bath does not depend on the particular choice of 
unitary that we perform, and furthermore the fidelity, 
being unitarily invariant, depends only on r\,„t. For the 
situation of a perfect cluster state (T = 0) with tg = 
and p;^^ = |+)(+| we have 

^(ft„,P„J = i(l + e-^"*cosAt), (13) 
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FIG. 1: Contour plot of fidelity as a function of coupling a, 
to a zero temperature bath and measurement time t. We set 
A = 1 and show ten equally spaced contours between F = 
(a = 0, t = tt) and F = 1 {a = Q,t = 2nn). Each shaded band 
corresponds to an interval of 10%, for example, the white re- 
gions correspond to a fidelity of 90 < F < 100%, centered 
around multiples of r, the uppermost large grey region corre- 
sponds to 40 < F < 50% while the black regions correspond 
to < F < 10%. 

which is plotted in Fig. [T] For a fixed a, the local max- 
ima for fidelity occur slightly before the times given by 
multiples of t = 27r/A, due to the presence of the expo- 
nential factor, however we note that the analysis derived 
from the master equation (|4]) will only be valid for weak 
bath couplings a ^ A. 

We see that to obtain high fidelity, we perform the 
measurements at times given by multiples of r. In the 
large t limit the output state decoheres to the maximally 
mixed state which reflects a return to the pure clus- 
ter state. We also note that the evolution of the encoded 
quantum information between time and time to is dis- 
tinct from the evolution between tg and to -f t and we 
will show that, in general, the latter form of evolution is 
the typical way in which fidelity is lost. For comparing 
results here with those obtained in the general decoher- 
ence situation, we note that the measurement on qubit 
1 produces a Hadamard transformation of the encoded 
state, and consequently swaps the x and z components 
of the Bloch vector. 



D. Optimal Cooling Rate 

In any experimental realization of MBQC on a strongly 
coupled system, there will be a residual thermal coupling 
to an ambient background (the environment), at temper- 
ature T^g. Typically, this environment is warm compared 
to the relevant energy scale in the system, i.e. fcTbg ^ A. 
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The coupling to this background can be reduced, for ex- 
ample by screening the system from thermal noise, but 
usually it can not be eliminated altogether. The purpose 
of the cooling bath (at temperature Tboth) is to counter- 
act this residual heating effect, by actively cooling the 
system such that the lattice of spins is prepared in a 
highly entangled cluster state at an effective temperature 
kT <^ A, fcTbg. However, as described in the previous 
section, the coupling to this bath also has an unwanted 
effect, which is to reduce the fidelity of MBQC on the sys- 
tem, by disrupting the state of the system over the course 
of the computation. A reasonable question to ask, then, 
is how the fidelity of a calculation varies as the strength 
of the coupling to the cooling bath is varied. 

The effects of a cooling bath plus high-temperature 
background may be modeled by including separate Lind- 
blad terms for each of the baths in the master equation: 

i 

+Y,[a,^M+).{-\]m+ix.M-)-.{MW)) 

i 

+Y,{c^^M\+h{-mt) + p.M\-h{+\]m) , (14) 

i 

where atath and /3bath describe the coupling to the cool- 
ing bath at temperature Tbath, and and /3bg are the 
corresponding coupling strengths to the background en- 
vironment at temperature Tbg . For simplicity, wc assume 
that the background is very warm compared to the en- 
ergy gap in the system, fcTbg ^ A, and use ^ to deduce 
that abath = /^bath = 7, and also that the cooling bath is at 
a very low temperature fcTbath ^ A, such that /3bath = 0. 
In this limit, the master equation becomes 

i i 

+ 7E(^[|+)^M]^(*) + n-H+wm) ■ (15) 

i 

(Note that the effect of a non-zero temperature cooling 
bath can also be described by this master equation by a 
suitable redefinition of ttbath, /3bath and 7). 

To understand the effect of cooling on the fidelity of 
a computation, we consider the three qubit example of 
Sec. Ill C[ using the master equation p5)) . We assume 
that the system is initially in equilibrium under psp and 
that the measurements are performed on qubits 1 and 2 
at times t = Q and t = 27r/A. Between the measurements 
the system evolves according to the master equation (|T5l) , 
after which we calculate the fidelity between the actual 
output state on qubit 3 and the ideal output state. The 
behaviour of the fidelity function of the coupling 
Q!bath, for various values of 7, is shown in Fig. [2l 

From Fig. [2]it can be seen that, for large Ofbath, there 
is a high loss in fidelity due to rapid dynamics in between 
the measurement steps, that try to bring the system back 
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FIG. 2: Fidelity as a function of the cooling constant Ofbath 
for couplings to the background (from the bottom to the top) 
7 = 10~^, 10~^, lO""*, and 10~^. For each there is an optimal 
cooling rate that maximizes the fidelity. 

to equilibrium. At the other extreme, for a weak coupling 
to the cooling bath ctbath 0, the initial state of the sys- 
tem is highly mixed, due to the coupling to the warm 
background, and so the fidelity is also reduced. There is 
consequently a trade-off in terms of cooling strength, be- 
tween counteracting the heating effects of the background 
and reducing errors due to dynamics between measure- 
ments for MBQC on the system. Thus given any ambient 
background there is an optimal coupling of the system to 
the cooling bath. Provided the coupling to this cooling 
bath is under the control of the experimentalist, the op- 
timal coupling should be selected in order to maximise 
the fidelity of computations. 

Note that the two Lindblad terms in (|15p can be ab- 
sorbed into a single term, such that the effect of the two 
baths is the same as coupling to a single bath with a = 
cKbath + 7 and /3 = 7, so that using ([6]), the effective tem- 
perature of the bath is given by kT = A/ log(l-|-Q;bath/7)- 
In the subsequent sections we treat the system as if it 
were coupled to a single bath parameterized by a and [3. 



III. GENERAL DECOHERENCE IN MBQC 

For our simple X-rotation gate on a 3 qubit state, we 
found that the loss in fidelity of the encoded qubit takes 
on a particularly simple form. In this section, we general- 
ize this result for an arbitrary sequence of measurements 
in a MBQC scheme, performed at the multiples of the 
natural timescalc r. Within a general MBQC scheme 
on a d-dimensional lattice, one dimension is identified as 
"time" and a (c? — l)-dimensional logical state evolves 
through the lattice via measurements. We show that the 
decoherence of this logical state as MBQC proceeds along 
the time direction coupled via Hj to a bath at a given 
temperature is described by a single quantum operation, 
acting on the logical state, producing anisotropic deco- 
herence towards the maximally mixed state. The impor- 
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tance of this result is that the error model for the logical 
qubit is Markovian when we restrict to measurements at 
multiples of r on the cluster state. This error model in 
turn allows for the application of standard fault-tolerant 
thresholds. 



A. One-dimensional lattice 

We begin by considering single-qubit unitaries through 
MBQC on one-dimensional lattices, and will consider the 
general case in the next section. On a line, with qubits 
labeled sequentially left to right, consider the situation of 
already having performed — 1 projective measurements 
after a time {N — l)r, where r = 2tt/A is the natural 
measurement time. Consequently, the qubit at site N — I 
is in a pure state, while the qubits i < N — 1 are partially 
entangled having evolved back towards the cluster state 
under the full Hamiltonian for the spin-lattice system 
coupled to the thermal bath. The qubits i > iV — 1 are 
in an entangled state and are still awaiting measurement. 

Site N-1 

• • • -• • • o o • 

We map the qubits N, N + 1, N + 2, ■ ■ ■ to a system of 
unentangled qubits by applying the unitary 

Vn=Y[ f/«-i,z , (16) 

i>N 

where Ui-i^ is the controlled- unitary between qubit 
i — 1 and qubit i. This map transforms the Hamiltonian 
i/tot ^ if'^j = VAriJtotV^ such that qubits i > N - 1 
are uncoupled, and more importantly localizes the logical 
state to qubit N, 

Disentangle qubits to the right 




• • • • • • o • • 

\ 

Logical state 

after N- 1 measurements 

The dynamics of the original state is determined by map- 
ping under Vn, evolving under H[^^ and then mapping 
back with V^^. However, because the qubits at sites i > N 
are assumed to be in equilibrium, they are static under 
the dynamics and may be ignored, and so we only have 
to consider the dynamics of the measured qubits together 
with the logical state at qubit N. The logical state, local- 
ized at site N, decoheres during the time interval {N—1)t 
to Nt under H' , where the qubit interacts with site A^ — 1 
through the cluster Hamiltonian term —^Zn-iXn, and 
with the bath through the XNjZjyqjyj. 

The state of the measured qubits and qubit N evolves 
according to the map £p'*{p) = J2i^i,TpWj ^, where 



the Kraus operators Wj for the entire system can be ex- 
pressed in terms of the Kraus operators Mi of Eq. ([7]) 

as 

W^/={»l,...,««},r = t/(Af,,,^®M,,,^0- • -^M.^^rW^ (17) 

and where U is the product of all the controUed-Z uni- 
taries for bonds to the left of site N. 

The evolution has the effect of partially entangling the 
logical state at N with the other qubits. The input state 
for the A^"" projective measurement is then 

Tri N-i[£rip)] = 1(1 + ■ S) , (18) 

and we obtain the components of its Bloch vector rjvr 
via "Tnt = Tr[(7£'^°'(p)]. Consequently, we can determine 
the equations of motion for this vector as a function of 
time from the master equation. 

For simplicity, we go to the interaction picture, with 
pi{i) = eyi\)[~iHct\p{t) eyiY>[iHct\ and obtain the equa- 
tions 

dt{XN{t))i = (« - /?)(Zw„i(i)), - ^{a + l3){XN{t))i 
dt{YN{i))i = -(a + P){YN{t))i 

dt{ZN{t))j^-]^{a + P){ZN{t))i, (19) 

where {A{t))i ~ Tv[pi{t)A\ for any observable A and a 
and (3 the temperature dependent coupling parameters 
as in (|4]). The equation for the Z component holds for 
any site s and so {Zs{t))i = exp[— ^(a -I- /3)t]. Now all 
computational measurements in MBQC on the cluster 
take place in the X-Y plane, and so {Zs{t)) = for all 
s < N and for all time t after site s has been measured. 

The components of the logical state in the interaction 
picture then evolve as 

3 

{XN{t))i = XNCxp[--{a + I3)t] , 
(iw(i))/ =2/Jvexp[-(a + /3)t], 

{ZN{t))i^ZNCxp[-^{a + m- (20) 

However, at times t = nr, we have pi{t) = p{t), and so 
for these times the decohcrence to the maximally mixed 
state is deduced from the interaction picture results, and 
agrees with the explicit example of the 3 qubit system in 

Sec, lire] 

The result of this analysis is that the MBQC scheme 
along the line of qubits with a free Hamiltonian and in 
contact with a bath at a finite temperature can be de- 
scribed in simple terms for any sequence of measurements 
on the logical state at times which are multiples of t, us- 
ing a fixed Markovian noise operator T. For a sequence 
of measurements • ■ • 0Af} in the X-Y plane la- 

belled at each site by an angle from the X axis and 
with outcomes {si, S2, . . . , sjv}, the single qubit logical 
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state is processed as 



(21) 



where for any A we denote [A][p] = ApA\ H = |+)(0| + 
|— X{(l)) = exp(— and Tip] is the quantum 
operation given by 



^[p] = Pip + P2XpX + p^YpY + piZpZ , 



(22) 



with 



pi = i(l + w)(l + u;2) p2 = 

P3 = i(l-w)(l + u;2) p4 = i(l + u;)(l_u.2) (23) 

and w = exp[— (a + f3)t/2]. 

It is also clear from this analysis why in the case of an 
arbitrary X rotation performed with three qubits, that 
the first evolution is slightly different from the second 
one: before the first measurement there are no qubits to 
the left of the first site to affect the logical state, while 
the qubits to the right are already in their equilibrium 
state, and so the evolution of the state is localized to the 
first site. 



B. General lattices 

The analysis of the last section can be extended to 
higher-dimensional lattices. For example, if upon local- 
ization to a site s using an analogous unitary to Vn, the 
single qubit logical state has k neighbouring sites, labeled 
1, 2, 3, . . . , fc, then generalizes to 



dt{X,it))i 

dt{Ys{t))i 
dt{Zs{t))i 
dt{Zi---Zk{t))i 



{a^ l3){Zi---Zk{t))i 
^ + ^a + /3)(X,(t)), 

{a + p){Y,{t)), 



2 

k + l 



{a + l3){Z,{t))i 

(a + /3)(Zi---Zfe(t))j 



(24) 



For general MBQC on for example the 2D or 3D lat- 
tice, we assume that the adaptive measurements are per- 
formed in steps, with a time interval r before the next 
round of measurements. Each round of measurements is 
composed of a set of Z measurements to eliminate qubits 
from the lattice and a set of measurements in the X-Y 
plane to propagate correlations and to perform the de- 
sired computational transformations. 

After TV — 1 such measurement steps, we may formally 
disentangle the qubits to be measured at steps N,N + 



Measured 
qubits 




Disentangled 
qubits 



FIG. 3: Schematic diagram for the localization of the logical 
state to the qubits Qm after A*' — 1 timesteps. The qubits 
to the right of Qn are disentangled with a unitary Vjv, while 
those to the left have already been measured. The qubits 
Qn-1 are in a pure product state, having just been measured. 



1, TV -I- 2, . . . and transform to a system where the logical 
state is localized to a set of qubits, which we denote Qn- 
Once again, the qubits to measured at stages iV -I- 1, -|- 
2, . . . are in their equilibrium states and can be ignored 
for the timestep. The reduced state for the logical state 
after a time r is determined from terms of the form (M)/, 
where M is a product of Pauli operators on Qn and its 
surrounding qubits. 

A site i on the lattice contributes to the equation of 
motion of a general observable M according to the follow- 



ing rule: it contributes 



{M)i if either [M, Z,] 



{M, Ki] = or if {M, Z,} = {M, K^} = 0; it contributes 
-{a + (3){M)I + {a-(^){MK^)I if{M,Z,} = [M,K^] = 0; 
and zero otherwise. Consequently, awkward terms can 
arise when Af contains an X observable. These terms 
couple the equations of motion for observables on Qn 
with observables on Qn plus its neighbours, however 
for MBQC on the cluster state these equations decou- 
ple, from the following argument. 

For an observable M on Qn containing m Xi observ- 
ables, its equation of motion will only have terms of the 
form {M)j and {MKi)j. However MKi is an observable 
with (to— 1) X's on Qn and a number of Z^s in Qn-i, the 
set of qubits that have just been measured. If we iterate 
and obtain the full set of coupled equations that deter- 
mine (MKi) J we arrive at a dependence on observables 
without any X observables and with at least one Zs on 
site s in Qn~i- The equation of motion for such an ob- 
servable M' is of the form dt{M')i = -p{a + (3)/2{M')i 
for some integer p. Furthermore, if s was measured in the 
X-Y plane then initially {M')j = and so will remain 
zero for the whole time interval. Retracing the chain of 
coupled equations we find that each problematic term of 
the form {MKi)j vanishes for t > and the equations 
of motion for the observables on Q^r are decoupled pro- 
vided each qubit in Q^r has at least one neighbour in 
Qn-1 that was measured in the X-Y plane. 

The expectation of an observable M on Qn will evolve 
as {M)i = Moe"''-""'"'^'*/^ for some integer q, and for 
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Pauli observables Mi, M2 . . . , Mj- on sites 1, 2, . . . , fc we 
have that \{M^ ■ ■ ■ Mi,{t)) j\ > |(Mi(<)),| • • • |(A4(0)/|, 
with equahty coming when the sites do not share any 
neighbours. 

C. Fault tolerance 

With a simple Markovian description of the errors 
present in our scheme, we can consider fault-tolerant 
MBQC. There are two sources of errors in the dynam- 
ical setting that we are considering. First, the equilib- 
rium state for the system is at a non-zero temperature, 
and so there are preparation errors due to an imperfectly 
prepared cluster state. Second, errors occur due to the 
dynamics between measurements, and can be viewed as 
storage errors on the qubits for a given timestep. For 
sufficiently low rates, MBQC on a 3D lattice has been 
shown to be fault-tolerant for both of these sources of er- 
rors [H, [IS] • If the state distillation protocol of Ref. [l^l 
is used, the error threshold is set by the bulk topological 
part of the error correction scheme, which in turn can 
be related to a phase transition in the classical random 
plaquette gauge model [2l|. 

Our initial state is a thermal state, static under the 
dynamics, prepared by cooling with the bath. Such a 
thermal cluster state at a temperature T is obtained by 
applying Z errors to a perfect cluster state with proba- 
bility = (1 + exp(A/(fcr)))-i. 

For the cubic lattice model, the dynamics in between 
measurement steps produce an error channel on the indi- 
vidual qubits no worse than a quantum operation of the 
same form as but with coefficients 

= i(l + u;)(l-f p2 = ^(l-«^)(l-^«') 

P3 = \{l - w){\ + p4 = ^(1 + w){l- w') , (25) 

and with w = e^'"'^'^^'^/^. Consequently, the resultant 
errors for a cubic lattice are no worse than those obtained 
by application of the local depolarizing channel T[p\ = 
(1 - ps)p + ^{XpX + YpY + ZpZ) with = + 
w)(l — w^), on each individual qubit. 

The combined effect of these two errors leads to inde- 
pendent errors on each qubit in the lattice with effective 
parameter q = Pp„p -I- |ps (c.f. [3). The threshold for 
such errors is given by g < 0.0293 [2l|. Thus, if errors 
due to the dynamics can be neglected, i.e. when ps 0, 
the error threshold for preparation errors corresponds to 
a temperature bound of T 0.28A. Conversely, if errors 
due to preparation can be neglected, Pprop —* 0, the error 
threshold corresponds to a threshold for the environmen- 
tal couplings of (a + P)/A ^ 4.6 x 10~^. 

If this environment consists of a infinite temperature 
background parametrized by 7 and a zero-temperature 
cooling bath parametrized by ofbath as in Sec. Ill Dl the pa- 
rameter q is a function of these two parameters. The con- 
straint g(abath,7) = 0.0293 defines the threshold value 



of 7 implicitly in terms of abath- We may then maxi- 
mize this 7 over the bath couplings and deduce an over- 
all threshold of 7/A ~ 3.4 x 10"'"' for the coupling to the 
environment provided that the cooling rate for the bath 
is set at a "Goldilocks value" of abath/A ~ 2.26 x 10^^. 
For this cooling rate the system is not so hot that large 
preparation errors destroy entanglement, and it is not 
so cold that large storage errors erase the logical state. 
That is, if the coupling to the background environment 7 
is below this threshold, it is possible to devise a cooling 
bath that allows for fault-tolerant MBQC. 



IV. DISCUSSION 

When the ground state of a physical system provides 
a resource state for MBQC we must necessarily take into 
account the system's dynamics. As we have discussed 
there are several sources of complication compared with 
MBQC on a static resource. While we may prepare the 
system very close to its ground state, any measurements 
we then perform on it will produce excitations and for a 
general adaptive measurement scheme, involving classical 
feed-forward, the resultant dynamics between measure- 
ments will perturb the state and affect the computation. 

For the simple, dispersionless Hamiltonian ([T|) describ- 
ing a lattice of spins we showed that measurements 
should be performed at a characteristic clock speed 2tt/A 
defined via the energy gap A. However, the presence of 
environmental interactions further complicates matters. 
For the environment, we considered both ambient back- 
ground effects and also the effects of a thermal bath used 
to prepare and maintain the lattice system. We found 
that an optimal cooling exists, which is a trade-off be- 
tween adequate shielding of the system from a hot back- 
ground and providing slow dynamics that allow adaptive 
measurements. Furthermore, the loss in fidelity due to 
this dynamics is conveniently described in terms of a sin- 
gle quantum operation that acts on the logical state. 

The importance of our results is that under certain con- 
ditions, the environment produces Markovian errors on 
the logical state and is thus amenable to error-correction. 
Our results are general and do not depend on the type 
of lattice or its dimensionality. In the particular case of 
a cubic lattice we may invoke fault-tolerance results for 
MBQC in the presence of local independent depolarizing 
errors to obtain a threshold of T ~ 0.28 A for the tem- 
perature of the prepared state when dynamics may be 
neglected, and a threshold of {a + P)/A ~ 4.6 x 10""^ 
for the ratio of environmental couplings to energy gap 
when the storage errors dominate. In addition, we ob- 
tained a threshold of 7/ A ~ 3.4 x 10~^ for the coupling 
to a high temperature environment provided there is a 
zero-temperature cooling bath with coupling abath/A ~ 
2.26 X 10"'^ to the lattice system. 

Several issues remain that deserve investigation. For 
example, we have not discussed possible imperfections in 
the Hamiltonian or measurement errors, both of which 
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would modify the above thresholds. Furthermore, the 
free Hamiltonian behaviour suggests the obvious strategy 
of performing all measurements at or near the clock cycles 
of r = 27r/A. However more complicated measurement 
strategies may exist that produce high fidelities in the 
presence of a fixed cooling. 

While the above formalism may be adapted to different 
settings or more particular questions, another key out- 
standing issue is the effect of finite-time measurements 
in which the measurements themselves are not instanta- 
neous but are spread over some small finite interval of 
time. Such a situation requires a more elaborate analysis 
than the one presented here, especially when the mea- 



surement time becomes comparable with the clock cycle 
time r. 
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